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ABSTRACT 

With recent Chandra observations, at least 75 percent of the X-ray background in the 2 — 10 keV energy 
range is now resolved into discrete sources. Here we present deep optical, near-infrared, submillimeter, 
and 20 cm (radio) images, as well as high-quality optical spectra, of a complete sample of 20 sources 
selected to lie above a 2 — 10 keV flux of 3.8 x 10~^^ erg cm~^ s~^ in a deep Chandra observation of the 
Hawaii Deep Survey Field SSA13. The 13 galaxies with / < 23.5 have redshifts in the range 0.1 to 2.6. 
Two are quasars, five show AGN signatures, and six are z < 1.5 luminous early galaxies whose spectra 
\ show no obvious optical AGN signatures. The seven spectroscopically unidentified / > 23.5 sources have 

^ ■ colors that are consistent with evolved early galaxies at 2; = 1.5 — 3. 

Only one hard X-ray source is significantly detected in an ultradeep submillimeter map; from the 
submillimeter to radio flux ratio we estimate a millimetric redshift in the range 1.2 — 2.4. None of the 
remaining 19 hard X-ray sources is individually detected in the submillimeter. These results probably 
reflect the fact that the 850 /im flux limits obtainable with SCUBA are quite close to the expected fluxes 
■ from obscured AGN. The ensemble of hard X-ray sources contribute about 10% of the extragalactic 

^ ' background light at submillimeter wavelengths. 

, From the submillimeter and radio data we obtain bolometric far-infrared luminosities. The hard X-ray 

sources have an average ratio of bolometric far-infrared to 2 — 10 keV luminosity of about 60, similar to 
that of local obscured AGN. The same ratio for a sample of submillimeter selected sources is in excess 
of 1100; this suggests that their far-infrared light is primarily produced by star formation. 

Our data show that luminous hard X-ray sources are common in bulge-dominated optically luminous 
\ galaxies with about 10% of the population showing activity at any given time. We use our measured 

. bolometric corrections with the 2 — 10 keV extragalactic background light to infer the growth of super- 

1-^ ' massive black holes. Even with a high radiative efficiency of accretion (e = 0.1), the black hole mass 

density required to account for the observed light is comparable to the local black hole mass density. 

O \ Subject headings: cosmology: observations — galaxies: distances and redshifts — galaxies: evolution — 

^ , galaxies: formation — galaxies: active — quasars: general 
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1. INTRODUCTION (Schmidt et al. 1998). Thus, an additional population of 

»r, ii, oc ( ■ i- ■ 1 j-i, • • either absorbed or flat spectrum sources is needed to make 

After more than 35 years of intensive work, the origin 111 i 1 • i 

r-i-i,ujv 1,1 J ivo'o\ ■ i-ii 4. r 11 up the background at higher energies. 

„^ of the hard X-ray background I XRB) IS still not fully un- ^^t^t^ ,1 • i 1 , . , -,1 • .1 n 

. : ■ , , 1 rri, VOID \, i- - 4- -i ni-c^x -4.1, -4. XRB synthesis models, constructed withm the frame- 

H derstood. Ihe XRB photon intensity, FIE), with units , ^ .V^^t ■• 1 i 1 i , 

d ' I, . -2 -1 1 \r-i -1 u ■ 4- J u work of AGN unifacation schemes, were developed to ac- 

. . 1 photons cm s kcV sr , can be approximated by a . r .^ , 1 ■ . •, r ,1 ^rm-. 1 ; 1 • 

1 Aip-r \, • 4.1, I, i- count for the spectral intensity of the XRB and to explain 

power-law, F(h) — Ah ^ , where £/ is the photon energy ; • ;i 1 , , n i , 

in keV. The HEAOl A-2 experiment (Marshall et al. 1980) ^"^^^ w u ?n^n !1 ^^1. '"^..'''^''''Fx. t 

found that the XRB spectrum from 3 - 15 keV is well de- f ^oltjer 1989; Madau, Ghise Imi & Fabian 

.11, 1 , • J 1 A J 4. 1,- 14. u u 1994; Matt & Fabian 1994; Comastri et al. 1995; Zdziarski 

scribed by a photon index 1 ~ 1.4, and this result has been , ' „. n i ;■ t^t-i 

n J J 4- J J 4- 1 • u 4- 1 et al. 1995; Gilh, Risaliti, & Salvati 1999; Wilman & 

conhrmed and extended to lower energies by recent analy- ^r- ■■ tt ■ on, ■, t 

ses of ASCA (Chen, Fabian, & Gendreau 1997; Gendreau Fabian 1999; Miyaji Hasmger & Schmid 2000) In 

et al. 1995; Miyaji et al. 1998; Ishisaki et al. 1998) and """"'^f orientation of a molecular torus 

BeppoSAX (Vecchi et al. 1999) data. At soft (0.5 - 2 keV) surroundmg-^the nucleus determines the classification of 

X-ray energies, 70 to 80 percent of the XRB is resolved ^ source. The models invoke along with a population 

• 4_ , u 4.U u/^oArr 4. 11-4. ^ ' 4. ot uuobscurcd type-1 AGN whose emission from the nu- 

mto discrete sources by the i^CD^i satellite (Hasmger et , ,. ,, , , ,• , , ;• r - . ■ 

1 1nno^ 4- £ 4-u 4-- 11 -J 4.-C J cleus we see directly, a substantial population of mtrm- 

al. 1998). Most of these sources are optically identifaed • ,, , , a ^t^t , , , , , 

, J 4.- 14.- 1 • /\r<-\r\ -4.1, 4. sically obscured AGN whose hydrogen column densities 

as unobscured active galactic nuclei (AGN) with spectra ^ .,_ni ^^nr , ; , , , , , 

4., 4. 4. 4. , 4-i-4.u a4- vDTD 4. of A// ~ 10 — 10^ Cm ^ arouud the nucleus block our 
that are too steep to account for the nat XRB spectrum 
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line-of-sight. 

A significant consequence of tlie obscured AGN mod- 
els is tliat large quantities of dust are necessary to cause 
the obscuration. The heating of the surrounding gas and 
dust by the nuclear emission from the AGN and the sub- 
sequent re-radiation of this energy into the rest-frame far- 
infrared (FIR) suggests that the obscured AGN should 
also contribute to the source counts and backgrounds in 
the FIR and submillimeter wavelength regimes (Almaini, 
Lawrence, & Boyle 1999; Gunn & Shanks 2000). 

Due to instrumental limitations, the resolution of the 
XRB into discrete sources at hard energies had to wait 
for the arcsecond imaging quality and high-energy sensi- 
tivity of Chandra. Deep Chandra imaging surveys are now 
detecting sources in the 2 — 10 keV range that account 
for 60 to 80 percent of the hard XRB (Mushotzky et al. 
2000, hereafter MCBA; Giacconi et al. 2001; Garmire et al. 
2001), depending on the XRB normalization. The mean 
X-ray spectrum of these sources is in good agreement with 
that of the XRB below 10 keV. Furthermore, because of 
the excellent < 1" X-ray positional accuracy of Chandra, 
counterparts to the X-ray sources in other wavebands can 
be securely identified (MCBA; Brandt et al. 2000; Giacconi 
et al. 2001). In this paper we determine spectroscopic red- 
shifts and optical, near- infrared (NIR), submillimeter, and 
radio properties of a complete hard X-ray sample drawn 
from the MCBA deep Chandra observations of the Hawaii 
Deep Survey Field SSA13. We take Ho = 65 /les km s"i 
Mpc~-^ and use a f^M — \^^t^ — \ cosmology throughout. 

2. SAMPLE AND OBSERVATIONS 

The present study is based on a 100.9 ks X-ray map 
of the SSA13 field that was observed with the ACIS- 
S instrument on the Chandra satellite in December 
1999 and presented in MCBA. The position RA(2000)= 
13'' 12" 21.40^ Dec(2000)= 42° 4l' 20.96" was placed 
at the aim point for the ACIS-S array (chip S3). Two 
energy-dependent images of the back-illuminated S3 chip 
and the front-illuminated S2 chip were generated in the 
hard (2 - 10 keV) and soft (0.5 - 2 keV) bands. MCBA 
chose the hard band energy range to be 2 — 10 keV to facil- 
itate comparisons with ASCA data. In the present paper 
we likewise use the 2 — 10 keV range. Other recent Chan- 
dra studies of the XRB have used either the 2 — 8 keV 
range (Brandt et al. 2000; Hornschemeier et al. 2000) or 
the 2 — 7 keV range (Fabian et al. 2000; Giacconi et al. 
2001) to minimize the backgrounds. 

We provide here a detailed description of the data re- 
duction techniques that were employed by K. Arnaud in 
December 1999 to analyze the SSA13 Chandra image for 
the MCBA paper. Improved X-ray data analysis tech- 
niques will be presented in Arnaud et al. (2001). 

The X-ray images were prepared with xselect and as- 
sociated ftools at GSFC. ACIS grades 0, 2, 3, 4, and 6 
were used, and columns at the boundaries of the readout 
nodes, where event select does not work properly, were re- 
jected. For the S3 chip the light curve was examined and 
times with high backgrounds were rejected, giving a total 
exposure of 95.9 ks; for the S2 chip the full time of 100.9 
ks was used. The images were examined in chip coordi- 
nates to identify and remove bad columns and pixels. For 
S3 the spectrum of the entire chip was extracted and the 



Si fiuorescence line was used to determine the gain. For 
the front-illuminated S2 chip radiation damage caused a 
systematic change in gain and spectral resolution across 
the chip. An observation from the Chandra Archive of the 
in-flight calibration sources was analyzed to determine an 
approximate correction for the spatial gain dependence. 
The PHA values were divided by (1-CHIPY*0.0002) to 
correct the gross gain variation. The remaining system- 
atic changes on the fluxes are small compared to statisti- 
cal uncertainties. The calibration sources were also used 
to determine the conversion from adjusted PHA to energy. 

Once the images in the 2 — 10 keV and 0.5 — 2 keV 
energy bands were generated for the chips, a simple cell 
detection algorithm was used to flnd the sources. The 
fleld area was stepped through in 2" steps, and the counts 
within a 5" diameter aperture (chosen to maximize the 
enclosed source counts throughout the image without be- 
coming substantially background dominated) were mea- 
sured, together with the background in a 5" — 7.5" radius 
annulus around each position. The average background 
was 4.7 counts in a 5" cell in the 2 — 10 keV S3 image and 
1.5 counts in the 2 — 10 keV S2 image. The distribution of 
counts is Poisson. A cut of 17 counts in the 2 — 10 keV S3 
image and 10 counts in the S2 image ensures that there is 
less than a 20% probability of a single spurious source de- 
tection in the entire sample. The background subtracted 
counts were searched for all positions within 4.5' of the op- 
tical axis that satisfied the appropriate count criteria. The 
source counts were next corrected for the enclosed energy 
fraction within the 5" aperture; the correction is small for 
this choice of off-axis radius and was determined from a 
second-order polynomial fit to the ratio of the 5" to 10" 
diameter aperture counts as measured from the brighter 
sources in the hard and soft band images. The final po- 
sitions were obtained with a centroiding algorithm that 
determined the center of light of the X-ray sources. 

The initial source positions were determined from the 
Chandra aspect information, but with a plate scale of 
0.4905" per pixel. The absolute pointing was then refined 
using the 10 sources with 5<t radio counterparts (see § 2.4). 
The offsets are 3.2" W and 0.3" N with no roll angle cor- 
rection. With the offsets applied, the dispersion between 
the radio and X-ray positions for the 10 sources is 0.4". 
The present positions should be more accurate than those 
given in MCBA. Figure |l| shows the Chandra hard band 
image of SSA13 with the source positions identified by the 
small circles. The large circle illustrates the 4.5' radius 
region used in this paper. 

The conversion from counts to flux depends on the shape 
of the source spectrum. For a power-law spectrum with a 
photon index F, the counts [photons s^^] in an energy band 

El to E2 are given by Ne^-e^ = « E'^dE, 
where A{E) is the effective detector area at energy E. 
Once the normalization of the spectrum, k, is determined 
from the observed counts, the flux in the 2 — 10 keV band 
is /hx = 1^° (FiE) E-^) dE where F = n E^-^ . The 
energy index is a = F — 1. 

Often F = 2 (typical of unabsorbed soft band sources) is 
assumed; however, in the absence of a correction for the ac- 
tual opacity, such a procedure underestimates the average 
conversion of counts to hard X-ray flux. An alternative 
approach is to determine the value of F for each source 
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from the ratio of soft to hard band counts, A^o. 5-2/^2-10- 
We adopt this approach here; however, because there is 
substantial uncertainty in the individual source F values, 
we use the counts- weighted mean photon indices of 1.2 
(hard) and 1.4 (soft) to determine the hard and soft band 
fluxes. Our procedure has the advantage of giving the cor- 
rect conversion for the ensemble of sources and therefore a 
correct comparison with the hard XRB, but it will result 
in errors in individual source determinations. The typical 
error is not large; for example, a source with an actual 
r = 2 spectrum would have its 2 — 10 keV flux overesti- 
mated by a factor of 1.35. A subsequent paper (Arnaud 
et al. 2001) will give the derived fluxes using the best fit 
spectral parameters to the individual sources. 

For the S3 chip, the final flux calibrations were made us- 
ing an array of effective areas versus energy at 12 positions. 
For the S2 chip a single conversion factor of 2.6 x 10~^^ 
erg cm~^ count ~^ was used in the 2 — 10 keV band and 
4.5 X 10~^^ erg cm~^ count"^ in the 0.5 — 2 keV band. 

To check our conversion of counts to flux, we ex- 
tracted the public ASCA and Chandra observations of 
G21.5, a compact crab-like supernova remnant with a well- 
determined simple spectrum. The vast difference in angu- 
lar resolution between ASCA and Chandra means that we 
can only make direct source flux comparisons for small 
{9 < 2') or point sources. We can also only make direct 
comparisons of data from instruments with different spec- 
tral resolutions for a simple spectrum source since there 
the transformation from source counts to flux is less model- 
dependent than it would be for a complex spectrum source. 
We derived the column density and power-law index for 
G21.5 from the Chandra S3 data and found that the values 
agreed well with the ASCA values. The measured Chandra 
flux was within 10% of the ASCA flux. At present there 
are no other time-stable compact or point sources with 
simple spectra in the public database that can be used for 
such a comparison. 

We selected all sources that were more than 15" from 
the chip edges and that have 2 — 10 keV fluxes greater 
than 3.8 x 10~^^ erg cm~^ s~^, corresponding to the above 
counts cut-offs and the appropriate calibrations. This de- 
tection threshold is slightly higher than the value quoted 
in MCBA (3.2 x 10~^^ erg cm^^ s^^) in order to obtain 
a complete flux-limited sample that is uniform over our 
4.5' radius area on the S2 and S3 chips. Table 1 details 
the 20 sources in the resulting 57 arcmin^ area, ordered by 
decreasing hard X-ray flux. The first six columns in Ta- 
ble 1 include the source identifications (the MCBA identi- 
fications are given in parentheses), RA(2000), Dec(2000), 
2-10 keV flux, 0.5 - 2 keV flux, and the value of F re- 
quired to match the ratio of the soft to hard X-ray fluxes 
in the absence of opacity. The remaining entries in Table 1 
are discussed in subsequent sections. 

While the low total X-ray source counts of most of 
our sources precludes a detailed analysis on an object- 
by-object basis, all of our sources are consistent with the 
Chandra PSF. They also all have rather hard X-ray spec- 
tra that are consistent with power-law or kT > 2 keV 
thermal spectra (Arnaud et al. 2001). Thus, the sources 
are not likely to be emission from groups of galaxies. The 
size constraints on the sources are not sufficiently restric- 
tive to discriminate between X-ray binaries, hot gaseous 
atmospheres, or AGN as the source of the X-ray emis- 



sion. We note that an improved comparison between the 
expected and observed images shows that object CXO 
J131159. 3-1-123928 (source 1 in the present sample), which 
was thought to be extended by MCBA, is in fact consistent 
with being a point source. 

2.1. Optical and Near-infrared Imaging 

We used the Low-Resolution Imaging Spectrometer 
(LRIS; Oke et al. 1995) on the Keck 10 m telescope in 
March 1997, February 1998 and 1999, and February and 
March 2000 to obtain i?-band and /-band images that 
cover the Chandra SSA13 field. 

Wide-field and deep HK' observations were obtained 
over a number of runs using the University of Hawaii Quick 
Infrared Camera (QUIRC; Hodapp et al. 1996) on the 
2.2 m University of Hawaii (UH) telescope and the 3.6 
m Canada-France-Hawaii Telescope. The if i^T' (1.9 ±0.4 
/im) filter is described in Wainscoat & Cowie (2001). The 
astrometry for the HK' image was established by fitting 
a linear solution to the 75 VLA sources (see § 2.4) with 
bright counterparts in the NIR image. The solution has 
a plate scale of 0.1891" per pixel and a rotation of —0.27 
deg from a standard N-E orientation. The dispersion of 
the radio and NIR positions is 0.43". Using a large over- 
lap sample of sources in the optical and NIR images, a 
third-order polynomial fit to the distortion in each LRIS / 
and B image was determined. Using the 90 VLA sources 
with bright / counterparts in the corrected image, we find 
a dispersion of 0.67" between the radio and optical posi- 
tions. 

The excellent < 1" X-ray positional accuracy permits 
the secure identification of the optical counterparts to the 
X-ray sources. Figure ^ shows thumbnail i?-band images 
of all 20 hard X-ray sources listed in Table 1. (Thumbnail 
/-band images can be found in MCBA.) In selecting the 
optical counterparts, we considered only / < 24.5 sources 
within a 1.5" radius of the nominal X-ray position. The 
optical separations are given in column 10 of Table 1. Six- 
teen of the sources have one such counterpart, and none 
has more than one. None of the offsets exceed 1", and 
the dispersion of the offsets is 0.5". Monte Carlo simu- 
lations with a randomized sample show that the average 
number of spurious identifications with an optical counter- 
part within 1" is 0.5. At 95 percent confidence, less than 
2 of the 16 identifications are spurious. For the 16 sources 
the magnitudes were measured at the optical center; for 
the remaining sources, the magnitudes were measured at 
the nominal X-ray position. For most of the sources the 
magnitudes were measured in 3" diameter apertures and 
corrected to approximate total magnitudes using an aver- 
age offset (Cowie et al. 1994); henceforth, we refer to these 
as corrected 3" diameter magnitudes. However, for the 
bright extended sources (7, 18, 19) we used 20" diameter 
aperture magnitudes and applied no correction; these mag- 
nitudes may be as much as a magnitude brighter than the 
corrected 3" diameter magnitudes given in MCBA. The 
B, I, and HK' magnitudes are given in columns 7, 8, and 
9 of Table 1. The la limits are approximately B — 27.6 
and / = 25.9. The la limits for the HK' magnitudes are 
not uniform over the field and are given individually in 
parentheses after the HK' magnitudes in Table 1. 

2.2. Keck Spectroscopy 
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We obtained high quahty optical spectra for 19 of the 20 
hard X-ray sources using LRIS sht-masks on the Keck 10 
m in March and April 2000. Source 14 was not observed 
because of mask design constraints. For the sources with 
I < 24.5 counterparts, we positioned the slit at the optical 
center. For the remaining sources, we positioned the slit 
at the X-ray centroid position. We used 1.4" wide slits 
and the 300 lines mm~^ grating blazed at 5000 A, which 
gives a wavelength resolution of ~ 16 A and a wavelength 
coverage of - 5000 A. The wavelength range for each ob- 
ject depends on the exact location of the slit in the mask 
but is generally between ~ 5000 and 10000 A. The obser- 
vations were 1.5 hr per slit mask, broken into three sets of 
0.5 hr exposures. Fainter objects were observed a number 
of times; the longest exposure was 6 hrs. Conditions were 
photometric with seeing - 0.6" - 0.7" FWHM. The ob- 
jects were stepped along the slit by 2" in each direction, 
and the sky backgrounds were removed using the median 
of the images to avoid the difficult and time-consuming 
problems of flat-fielding LRIS data. Details of the spec- 
troscopic reduction procedures can be found in Cowie et 
al. (1996). 

We successfully obtained redshift identifications for all 
13 sources brighter than / = 23.5 mag; the spectra are 
shown in Fig. 0, and the redshifts are given in column 11 
of Table 1. We classify the spectra into three general cate- 
gories: (i) quasars (broad-line sources), (ii) AGN (narrow 
and weak-line sources), and (iii) optically 'normal' galax- 
ies (no AGN signatures in the optical). Henceforth, we 
denote these categories by g, a, and n, respectively. We 
also denote spectroscopically unidentified sourc es by u and 
our one source with a millimetric redshift (see § |6.l[ ) by m. 

Two sources (sources 3 and 6) are the quasars previously 
known to be in the field (Windhorst et al. 1995; Campos 
et al. 1999). Their spectra are very similar, and both coin- 
cidentally lie at z = 2.565. These quasars are radio quiet. 

Five sources (1, 2, 5, 11, and 15) show emission line char- 
acteristics that may be indicative of AGN activity. Sources 
1 and 5 show [O II] , Ne III] , and weak Ne V] emission, along 
with Ca H and K and G-band absorption. Sources 11 and 
15 show narrow Lya and CIV emission. Source 2 shows a 
P-Cygni profile in Mg II and broad absorption in Fe II. 

Six sources (4, 7, 10, 12, 18, and 19) show no indication 
of an active nucleus in their optical spectra. We call these 
'normal' galaxies since they have absorption and emission 
line properties which are common in optically selected field 
samples. Source 4 shows Ha and weak [OH] and [QUI] 
emission and H/3 absorption. Source 7 shows Ha, [Oil], 
and [O III] emission and H/3 absorption. Source 10 shows 
narrow [Oil] emission and Mgll absorption. Source 12 
has weak [Oil] and H/? emission and strong [OIII] ab- 
sorption. There may be hints of Nelll] and NeV] . Source 
18 is rather unusual in that it has no Ha while Nil and 
S II are in emission. A high N II/Ha ratio has been used to 
classify objects as AGN (Keel et al. 1985), but the absence 
of Ha in source 18 is difficult to understand: in Veilleux 
& Osterbrock (1987) the highest ratio of Nil to Ha is 3:1, 
and no photoionization models (Ferland & Netzer 1983) 
have ratios larger than 2:1. Source 19 has Ha emission 
but otherwise only absorption features. 

Optically 'normal' X-ray luminous galaxies were 
thought to be relatively rare, unusual objects, perhaps ex- 
plained by beaming (Elvis et al. 1981; Moran et al. 1996; 



Tananbaum et al. 1997). They are hard to find by as- 
sociation since small X-ray error boxes are required to be 
certain of their identification (e.g., discussion in Schmidt et 
al. 1998). The very large surface density of such sources in 
our sample, ~ 400 deg^^, indicates that they are common. 
In fact, in our sample they are much more common than 
quasars. There are two plausible explanations for the lack 
of observed optical AGN characteristics: i) absorption due 
to dust and gas or ii) an actual lack of ultraviolet/optical 
emission, as is the case in many low luminosity objects 
(Ho et al. 1999). The line of sight column densities in- 
ferred from the X-ray spectra in § 3.1 are sufficiently large 
to obscure the optical AGN signatures, but the extent will 
depend on the geometry. 

2.3. Suhmillimeter Observations 

The suhmillimeter observations were made with the 
SCUBA instrument (Holland et al. 1999) on the James 
Clerk Maxwell Telescope. SCUBA jiggle map observations 
were taken in mostly excellent observing conditions during 
runs in February 1999 (7 observing shifts), February 2000 
(0.5 shift), and May- June 2000 (3.5 shifts). The maps were 
dithered to prevent any regions of the sky from repeatedly 
falling on bad bolometers. The chop throw was fixed at a 
position angle of 90 deg so that the negative beams would 
appear 45" on either side east-west of the positive beam. 
Regular "skydips" (Lightfoot et al. 1998) were obtained to 
measure the zenith atmospheric opacities at 450 and 850 
/im, and the 225 GHz sky opacity was monitored at all 
times to check for sky stability. The median 850 /im opti- 
cal depth for all nights together was 0.185. Pointing checks 
were performed every hour during the observations on the 
blazars 1308-f 326 or cit6. The data were calibrated using 
30" diameter aperture measurements of the positive beam 
in beam maps of the primary calibration source Mars and 
the secondary calibration sources OH231.8, IRC-l-10216, 
and 16293-2422. 

The data were reduced in a standard and consistent 
way using the dedicated SCUBA User Reduction Facility 
(SURF; Jenness & Lightfoot 1998). Due to the variation 
in the density of bolometer samples across the maps, there 
is a rapid increase in the noise levels at the very edges. The 
low exposure edges were clipped from our images. 

The SURF reduction routines arbitrarily normalize all 
the data maps in a reduction sequence to the central pixel 
of the first map; thus, the noise levels in a combined im- 
age are determined relative to the quality of the central 
pixel in the first map. In order to determine the absolute 
noise levels of our maps, we first eliminated the ^ 3cr real 
sources in each field by subtracting an appropriately nor- 
malized version of the beam profile. We then iteratively 
adjusted the noise normalization until the dispersion of the 
signal-to-noise ratio measured at random positions became 
~ 1. Our noise estimate includes both fainter sources and 
correlated noise. 

We centered on the positions of the hard X-ray sources 
and measured the submillimeter fiuxes using beam- 
weighted extraction routines that include both the posi- 
tive and negative portions of the chopped images, thereby 
increasing the effective exposure times. The 850 /^m sub- 
millimeter fluxes and la uncertainties are summarized in 
column 12 of Table 1. For most of the objects the la level 
is in the 1 — 2.5 mJy range; however, for the two sources 
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in the region where there is an uhradeep SCUBA image 
(Barger et al. 1998) the Icr detection threshold is 0.6 — 0.7 
mJy. 

2.4. Radio Observations 

A very deep 1.4 GHz VLA radio map of the SSA13 re- 
gion was obtained by Richards et al. (2001) using an 100 
hr exposure in the A-array configuration. The primary 
image covers a 40' diameter region with an effective res- 
olution of 1.6" and a Str limit of 25 micro- Jansky (/xJy). 
The radio fluxes were measured in 2.4" boxes centered on 
the X-ray positions; these fluxes are given in column 13 of 
Table 1. Of the 20 hard X-ray sources, 16 are detected in 
the radio above a Scr threshold of 15 /iJy, including 10 of 
the 13 sources with spectroscopic redshifts. The radio- X- 
ray offsets are given in column 14 for the 10 sources with 
5(7 radio detections within 1.5" of the X-ray source. The 
absolute radio positions are known to 0.1 — 0.2" rms. The 
dispersion between the radio and X-ray positions is 0.4", 
and the maximum separation is 0.7". The radio to opti- 
cal ratios for the hard X-ray sources are consistent with 
the sources being radio quiet. The radio properties will be 
discussed in more detail in Richards et al. (2001). 

3. X-RAY PROPERTIES OF THE HARD X-RAY SAMPLE 
3.1. Optical Depths 

With the exception of two 'normal' galaxies with F < 
— 1.5, the photon indices given in Table 1 range from —0.4 
to 1.82. The two quasars have T — 1.75 and T = 1.80, con- 
sistent with most unabsorbed soft band sources, whereas 
most of the remaining sources have photon indices that 
suggest substantial line-of-sight optical depths. 

If we generate counts- weighted mean photon indices for 
each population separately, we find 0.8 for the 'normal' 
galaxies, 0.9 for the AGN, 1.8 for the two quasars, and 
1.2 for the spectroscopically unidentified sources. The 
'normal' galaxies are presumably those where the AGN 
are the most highly obscured. Since the effective column 
density, N^ff, for observed- frame absorption in a source 
is related to the true hydrogen column density, Nh, by 
Neff ~ Nh/{1 + z)'^'^, flux corrections for absorption ef- 
fects become less important with increasing redshift. The 
spectroscopically unidentified sources are lik ely h igher red- 
shift analogs of the 'normal' galaxies (see § 4.2) but have 
softer spectra because of this redshift dependence of the 
absorption. 

In Fig. ^ we plot the logarithm of the ratio of the soft to 
hard X-ray fluxes versus redshift. Throughout the paper 
we use the notation of filled diamonds for quasars, filled 
triangles for galaxies with AGN signatures, filled squares 
for galaxies with apparently 'normal' spectra, and open 
squares for spectroscopically unidentified sources. We 
overlay on the data fixed Nh curves which we generated 
assuming an intrinsic F = 2 power-law spectrum and pho- 
toelectic cross sections computed for solar abundances by 
Morrison & McCammon (1983). Over the energy range 
0.5 — 7 keV, which determines the correction for these col- 
umn densities and redshifts, the cross-section, cr{E), can 
be well-approximated by a single power-law 

a{E) = 2.4 X 10^22 E-'^-^ cm"^ 
with E in keV. For all but one of the spectroscopically 



identified galaxies the X-ray flux ratios can be described 
by a rather narrow range of neutral hydrogen column den- 
sities from Nh = 2 x 10^^ cm~^ to 3 x 10^^ cm~^, although 
more sophisticated models with scattering could permit 
higher opacities. The Nh values and the true power-law 
indices are best determined directly from the X-ray spec- 
tra, as shall be discussed in a subsequent paper (Arnaud 
et al. 2001). 

Could the absence of sources with column densities 
above Nh = 3 x lO^'^ cm~^ be a selection effect? This 
is possible at low redshifts since the ratio of the absorbed 
to the actual 2 — 10 keV flux drops rapidly for N^ff > 10^^ 
cm~^; however, at high redshifts this N^ff corresponds to 
larger values of Nh than are observed (see Fig. ||). Even 
if we place all of the spectroscopically unidentified sources 
at z » 1, at most 3 of the 14 objects with z > 1 have 
column densities above Nh = 3 x 10^'^ cm~^. The sim- 
plest interpretation of Fig. |4| is that we are seeing most 
of the obscured AGN in the present sample. However, 
Compton-thick sources might be missed completely from 
the 2—10 keV sample, and these sources may be needed 
to explain the 30 keV peak in the XRB. Unfortunately, 
this issue cannot be decided until either more precise in- 
formation is obtained on the energy distribution of the 
individual sources, which, in principle, could deviate from 
power-law behavior, or until we can analyze the hardness 
of even fainter 2—10 keV sources. 



3.2. Redshift Distribution 

The redshift distribution of the spectroscopically iden- 
tified sample is shown in Fig. pi where we we plot the sur- 
face density of the sources in Kz = 0.5 redshift bins (filled 
squares) with Icr uncertainties. Below 2 = 3 the data are 
consistent with a constant surface density. If we also in- 
clude the unidentified sources spread uniformly through 
the z = 1.5to2; = 3 redshift interval (open squares), then 
the constancy with redshift becomes even more evident. 
The surface density is ^ 200 sources per square degree per 
bin or ^ 400 sources per square degree per unit redshift. 

The present data sample is too small to justify a de- 
tailed examination of the luminosity function. However, 
the redshift distribution of the hard selected sample is very 
similar to that of previous soft selected samples with sim- 
ilar limiting sensitivities. We illustrate this in Fig. || by 
overplotting the redshift distribution of the Lockman Hole 
i? 05"^ r sample of Schmidt et al. (1998) (filled diamonds). 
The limiting flux of this sample is 5.5 x 10~^^ erg cm~^ s~^ 
(0.5 — 2 keV) which, for a source with F = 2, would corre- 

^ ' 10 

-2 
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spond to a limiting flux of 4.7 x 10"^"' erg cm" 
keV), similar to our limiting flux of 3.8 x 10"^*^ erg cm 
s~^ (2 — 10 keV). The absolute surface density of the soft 
sample is far below that of the hard sample, indicating that 
most of the hard X-ray sources are substantially obscured; 
however, the redshift distributions are rather similar. The 
luminosity function evolution of soft X-ray samples has 
been extensively analyzed and shows a rapid rise between 
z = and z = 1.5, followed by relative constancy at higher 
redshifts (Miyaji et al. 2000). While a proper analysis of 
the evolution of the 2 — 10 keV luminosity function must 
await larger samples and a better understanding of the 
optical depths and K-corrections, the flrst impression is 
that the behavior will be similar to that inferred from the 
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Fig. 4. — Ratio of the soft to hard X-ray fluxes versus redshift. Quasars are denoted by filled diamonds, AGN by filled triangles, 'normal' 
galaxies by filled squares, and spectroscopically unidentifkal-sources by open squares. The latter, excluding source 8, are placed at 2 = 2; 
source 8 is placed at its millimetric redshift of 1.8 (see § p lf . The two quasars have been slightly offset in redshift for clarity, as have the 
unidentified sources 16 and 17. As illustrated in the figure, the observed fiux ratios can be described by an input F = 2 spectrum and a 
narrow range of Njj; curves are for Njj = 2 X 10-^^ cm~^ (soUd), Njj = 10^^ cm~^ (dotted), and Njj = 3 X 10^'^ cm~^ (dashed). 
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Fig. 5. — Filled squares with Icr uncertainties show the surface density of the spectroscopically identified hard X-ray sources in Az = 0.5 
redshift bins. For 2 = 3 — 4 we show the Icr upper limit as the downward pointing arrow. Open squares include the unidentified sources 
spread uniformly from z = 1.5 to 2 = 3. Filled diamonds show the surface density of the deep ROSAT soit X-ray survey of the Lockman Hole 
(Schmidt et al. 1998). 



soft samples but with a much higher normahzation of the 
luminosity function (about 7 times higher). 

3.3. Luminosities 

The intrinsic flux, Fint, is related to the observed flux, 

Fobs, by 

Pint = Fobs (1 + z)"^^^ 

For an unabsorbed spectrum with F = 2 the K-correction 
vanishes. An unabsorbed spectrum may be appropriate at 
the higher energies where opacity effects are not impor- 
tant. However, we believe that we can obtain a slightly 
improved estimate for the 2 — 10 keV luminosities by al- 
lowing for the average effects of the opacity as follows. In 
calculating our hard X-ray luminosities, we normalized the 
flux at 4 keV for the F = 2 spectrum to the flux at 4 keV 
calculated over the 2 — 10 keV energy range for a spectrum 



with counts- weighted mean photon index F = 1.2 (see § 
Then 

Lhx - 47r 4 (0.85) fnx 

where fnx is the 2 — 10 keV flux of Table 1 computed with 
the same F = 1.2 assumption. Our hard X-ray luminosi- 
ties are given in column 3 of Table If we had instead 
used F = 2, the computed 2 — 10 keV fluxes would be 
lower by a factor of 1.35 and our X-ray luminosities would 
be lower by a factor of 1.15. The X-ray luminosities of 
Fabian et al. (2000) and Hornschemeier et al. (2000) were 
based on F = 2; in later comparisons with their luminosi- 
ties, we ignore this difference since it is small compared to 
other uncertainties. 

Figure ^ shows our hard X-ray luminosities versus red- 
shift. Here the solid curve represents the detection limit. 
The open pentagons denote the ultraluminous starburst 
galaxy Arp 220, the ultraluminous, highly obscured AGN 
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Fig. 6. — Hard X-ray luminosities versus redshift for the thirteen I < 23.5 hard X-ray sources with spectroscopic identifications. Quasars 
are denoted by filled diamonds, AGN by filled triangles, and 'normal' galaxies by filled squares. The solid curve shows the detection limit. 
The open pentagons represent Arp 220 (least luminous), NGC 6240, and PC 1543+489 (most luminous). 



NGC 6240, and the radio quiet quasar PG 1543-f489, 
in ascending Lhx order (see § ^ for comparison. The 
hard X-ray luminosities of the spectroscopicaUy identified 
sources range from just over 10"*^ erg s^^ to ~ 10'*^ erg 
s~^. Sources at low redshift (z < 1) do not have the high 
X-ray luminosities of the sources at high redshift. To de- 
termine whether there are also low luminosity sources at 
high redshift will require Chandra observations that probe 
to much deeper flux levels. 

Even though the 'normal' galaxies, which typically fall 
near the detection threshold, are systematically less lumi- 
nous than the AGN and quasars, they are still extremely 
X-ray luminous; only the two lowest luminosity sources 
even overlap the local galaxy populations (Fabbiano 1989). 
Furthermore, the 'normal' galaxies are generally more lu- 
minous than the hard X-ray sources in the nuclei of nearby 
giant elliptical galaxies in the Virgo and Fornax clusters 
that were chosen for study on the basis of their black hole 
properties; these have 1 — 10 keV luminosities of 2 x 10""^ 
erg to 2 X 10*^ erg s"^ (Allen, Di Matteo, & Fabian 
2000). 

3.4. Contribution to the Hard XRB 



The integrated 2-10 keV light of the 20 hard X-ray 
sources in the sample corresponds to an extragalactic back- 



ground light (EBL) of 1.34 X 10 erg cm ^ s 



^ deg 



between 58 to 84 percent of the 2 — 10 keV XRB, depend- 
ing on whether we assume the HEAOl value of 1.6 x 10^^^ 
erg cm~^ s~^ deg~^ (Marshall et al. 1980) or more recent 
higher estimates of 2.3 x 10"^-^ erg cm~^ s~^ deg~^ (e.g., 
Vecchi et al. 1999). Throughout the rest of the paper we 
will conservatively adopt the Vecchi et al. (1999) value. 

If we include 6 of the 7 spectroscopicaUy unidentified 
sources in the z < 2 population (source 13 is the only 
source not seen in the _B-band and hence is the only very 
high redshift {z 3> 4) candidate), then 56% of the light 
arises from z < 2 sources. If we instead restrict to the 
spectroscopic sample, then 34% of the light arises from 



the z < 2 population. The four sources with spectroscopic 
redshifts z > 2 contribute 13% to the light, and source 13 
contributes just 1%. 

4. OPTICAL AND NEAR-INFRARED PROPERTIES OF THE 
HARD X-RAY SAMPLE 

4.1. Magnitudes 

At the time of this writing, all but one (source 13) of 
the hard X-ray sources was detected above the 2a level in 
the NIR image, all but two (13, 16) in the /-band image, 
and all but one (13) in the B-band image. Source 13 has 
since been detected in the NIR using the CISCO infrared 
camera on the Subaru 8 m telescope (Cowie et al. 2001). 

Figure |^ shows the redshift versus / magnitude distribu- 
tion of our spectroscopicaUy identified / < 23.5 hard X-ray 
sample (large symbols) and, for comparison, an optically 
selected / < 24 field galaxy sample (small symbols). Hard 
X-ray sources with 5(t radio detections are indicated by 
surrounding open boxes. The tracks are from Coleman, 
Wu, & Weedman (1980) for an early-type galaxy (solid), 
an early spiral galaxy (dashed), and an irregular galaxy 
(dotted) with absolute magnitudes Mj = —22.5 in the as- 
sumed cosmology. The 'normal' galaxies and AGN follow 
the upper envelope of the star forming field galaxy popu- 
lation. Thus, the hard X-ray sources predominantly lie in 
the most optically luminous galaxies. 

4.2. Colors 

Figures Ha, b show / — HK' versus redshift and / — HK' 
versus HK^ for 19 of the 20 hard X-ray sources (source 13 
is excluded). The overlays are Coleman, Wu, & Weedman 
1980 tracks for an early-type galaxy (solid curve) and an 
early spiral galaxy (dashed curve) with Mhk' = —25.0. 
The colors of the spectroscopicaUy identified z < 2 galaxies 
are in the range of the galaxy tracks; this is also consistent 
with their morphological appearance (see Fig. ||). Of the 
sources that were too optically faint for spectroscopic iden- 
tification, the NIR magnitudes and colors for four (sources 
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Table 2 



Number z Lax Lqpt L|.t^™ Lbol 

(lO"*" ergs s-i) (10*"^ ergs s^^) (10^'' ergs s'^) (lO''^ ergs s"^) (10*'' ergs s"*) 






2.000" 


1000 


480 


4500 


< 8800 


9400 


1 


1.048^" 


150 


41 


3600 


< 19000 


4300 


2 


1.320^^ 


230 


2400 


< 1200 


< 17000 


< 4600 


3 


2.565'^ 


930 


2.8 X lO'^ 


28000 


< 51000 


3.1 X 10^ 


4 


0.212'' 


1.9 


4.6 


< 18 


< 1800 


< 31 


5 


0.696^* 


28 


140 


800 


< 7900 


1000 


6 


2.565'* 


570 


49000 


11000 


< 41000 


62000 


7 


0.241" 


1.3 


8.0 


130 


< 2100 


140 


8 


1.800°^ 


120 


320 


6500 


5200 


7400 
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2.000" 


130 


270 


4700 


< 14000 


5500 


10 


1.427" 


59 


520 


1800 


< 12000 


2500 


11 


2.415^ 
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4400 


8200 


< 16000 


14000 


12 


0.585" 


6.5 


160 


220 


< 5100 


400 


13 


2.000" 


130 


94 


< 3200 


< 14000 


< 3800 


14 


2.000" 


120 


140 


6600 


< 12000 


7300 


15 


2.625=^ 


220 


2200 


< 5900 


< 10000 


< 9000 


16 


2.000" 


110 


110 


3600 


< 5600 


4200 


17 


2.000" 


110 


650 


16000 


< 9600 


18000 


18 


0.110" 


0.12 


23 


27 


< 680 


51 


19 


0.180" 


0.34 


14 


77 


< 1300 


93 




Fig. 7. — Redshift versus / magnitudes for the thirteen I < 23.5 hard X-ray sources with spectroscopic identifications (large symbols) and 
for the / < 24 field galaxies from the Hubble Deep Field and the Hawaii Deep Survey Fields SSA13 and SSA22 (small symbols). Quasars are 
denoted by filled diamonds, AGN by filled triangles, and optically 'normal' galaxies by filled squares. Sources with 5(t radio counterparts are 
enclosed in open squares. The hard X-ray sources are almost exclusively galaxies lying at or above L,, as indicated by the Coleman, Wu, & 
Weedman (1980) tracks for an early-type galaxy (solid), an early spiral galaxy (dashed), and an irregular galaxy (dotted) with Mj = —22.5. 



0, 9, 14, 16) suggest that they are early galaxies in the 
z > 1.5 redshift range (Fig. Crawford et al. (2000) 

also found a number of sources of this type in a sample of 
Chandra hard X-ray sources. We were able to use radio 
and submillimeter detections for one of the optically faint 
sources (source 8) to estimate a millimetric redshift with 
central value z = 1.8 (see § ^). We used radio detections 
and submillimeter la limits on sources 0, 9, 14, 16, and 
17 to estimate millimetric redshift upper limits of z —2.5, 



3.0, 2.4, 2.2, and 1.4, respectively. Sources 8 and 17 are 
somewhat bluer than the curves in Figs. ||, but their AGN 
might be contributing substantially to the rest-frame op- 
tical light. 

4.3. Bolometric Ultraviolet/ Optical Luminosities 

We estimate the AGN luminosities in the ultraviolet 
(UV)/optical by adopting a shape appropriate to the ra- 
dio quiet AGN (e.g., Zheng et al. 1997). At frequencies 
below the rest-frame Lyman limit we take the spectrum 
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Fig. 8. — (a) I — HK' colors versus redshift and (b) I — HK' colors versus HK' magnitudes for the 19 hard X-ray sources with optical 
and/or NIR detections (source 13 is not shown). Quasars are denoted by filled diamonds, AGN by filled triangles, 'normal' galaxies by filled 
squares, and spectroscopically unidentified sources by open squares. The latter, excluding source 8, are placed at 2: = 2; source 8 is placed 
at its millimetric redshift of 1.8. The overlays are the Coleman, Wu, & Weedman (1980) tracks for an early-type galaxy (solid curve) and an 
early spiral galaxy (dashed curve) with M^j^i = —25.0. The spectroscopically unidentified sources are likely early galaxies at 2 > 1.5 because 
of their red colors. 



to be a —0.8 power-law; this steepens to a —1.7 power-law 
at higher frequencies. Normalizing to the observed flux 
at the wavelength corresponding to the rest wavelength 
2500 A, /25oo(i+z)j the UV/optical luminosity is then 

LoPT = 47r dl (9.4 X lO^^) /25oo(i+.) (1 + •2)"' 

where is the luminosity distance in cm and /25oo(i+2) is 
in units erg cm~^ s~^ Hz~^. We have estimated /25oo(i-i-z) 
for the hard X-ray sample by interpolating from the ob- 
served fluxes corresponding to the magnitudes given in 
Table 1 and, where available, the U' magnitudes at 3400 
A from MCBA. For the more extended objects we have 
used the corrected 3" diameter magnitudes from MCBA 
rather than the 20" diameter magnitudes of Table 1 to ob- 
tain better limits on the AGN fluxes. Many of the sources 
may still have substantial galaxy light contamination, so 
LoPT is strictly an upper limit on the AGN contribution. 
The inferred bolometric UV/optical luminosities are given 
in column 4 of Table ^ and are plotted versus redshift in 
Fig. |. 

The four sources with spectroscopic redshifts beyond 
z — 2 are by far the most optically luminous sources in 
the sample. These are the two quasars (sources 3 and 6) 
and the two AGN (sources 11 and 15) with narrow Lya 
and CIV lines. The AGN are likely dominating the light in 
these sources at the observed optical and NIR wavelengths, 
as the observed colors are quite blue (see Fig. ^). 

5. AN OPTICALLY SELECTED SAMPLE: HARD X-RAY 
PROPERTIES AND CONTRIBUTION TO THE HARD XRB 

The foregoing section presented the optical nature of the 
hard X-ray sources. We now invert the approach and ask 
what are the hard X-ray properties of an optically selected 
sample. In particular, we would like to know what frac- 
tion of optical sources are significant contributors to the 
hard XRB and whether these are drawn from a particular 
subsample of the optical population. 



To address this issue, we use the complete subsample of 
1151 / < 24 galaxies and stars within a 4' radius of the 
optical axis of the Chandra pointing. The smaller radius 
was chosen to minimize aperture corrections to the X-ray 
counts so that the X-ray limits would be more uniform. 
For each of these sources, we extracted the X-ray counts 
in a 2" radius aperture centered on the nominal optical po- 
sition and converted them to 2 — 10 keV fluxes using the 
procedure outlined in § ^ The X-ray fluxes determined in 
this way are plotted versus / magnitude in Fig. 

Above a hard X-ray flux limit of 2.5 x 10~^^ erg cm~^ 
s~^, we recover the twelve / < 24 sources of Table 1 that 
are within the 4' radius. Monte Carlo simulations using 
an equal number of sources at random positions show that 
with this extraction aperture and flux limit, there will 
be, on average, 1.4 spurious cross- identifications of optical 
sources with X-ray sources (< 3 at the 95 percent con- 
fidence limit). Allowing for this contamination, we find 
that a fraction 0.009^Q gQ3 of the optical sources are X-ray 
sources, where the uncertainties correspond to the 68 per- 
cent confidence ranges. At a lower hard X-ray flux limit of 
1.0 X 10~^^ erg cm~^ s~^, we find 32 sources and expect 
an average contamination of 10.3 (< 16 at the 95 percent 
confidence limit). Here, a fraction 0.019 ± 0.005 of the 
optical sources are X-ray sources. 

Excluding the incompletely covered edges of the S3 and 
S2 chips, the observed area is 47 arcmin^, and the con- 
tribution to the 2 — 10 keV background of the 12 sources 
above the 2.5 x 10"^^ erg cm~^ s~^ hard X-ray flux limit 
is 9.5 X 10-12 grg cm-2 g-i deg-^ (41% of the hard XRB). 

If we assume that half of the sources in the flux range 
1.0 X 10--^^ erg cm-2 s-^ to 2.5 x 10"-^^ erg cm^^ s^^ are 
real, we estimate a further contribution of only 1.1 x 10-""^^ 
erg cm^2 s^^ deg^^; thus, increasing the X-ray sensitiv- 
ity adds little to the hard XRB contribution. Indeed, 
summing over all the optically selected sources above 
1.0 X 10-1^ erg cm-2 s'^ yields a total of 10.4±0.8 x IQ-^^ 
erg cm-2 deg-^, where the uncertainties are the 68 
percent confidence ranges based on randomized samples. 



10 




Redshift 



Fig. 9. — UV/optical luminosities versus redshift. Quasars are denoted by filled diamonds, AGN by filled triangles, 'normal' galaxies by 
filled squares, and spectroscopically unidentified sources by open squares. The latter, excluding source 8, are placed at z = 2; source 8 is 
placed at its millimetric redshift of 1.8. Source 16 has been slightly offset in redshift for clarity. 
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Fig. 10. — (a) Hard X-ray fluxes versus I magnitudes for the 1151 / < 24 galaxies and stars within a 4' radius of the optical axis. 
Above a hard X-ray flux limit of 2.5 X 10~^^ erg cm~^ s~^, the twelve / < 24 sources of Table 1 that are within this radius are recovered, 
(b) Hard X-ray fluxes versus absolute /-band magnitudes for the 172 7 < 23 spectroscopically identified sources within the 4' radius. The 
K-corrections have been computed using the Coleman, Wu, & Weedman (1980) SED for an early spiral galaxy. Five of the 25 galaxies with 
—22.5 < Mi < —24 have fnx > 1-0 x 10~^* erg cm~^ s~^. None of the 29 galaxies with Mj fainter than —20 is an X-ray source at this 
level. 



We conclude that most of the hard XRB from the op- 
ticaUy selected sources is dominated by the small number 
of sources with hard X-ray fluxes above 2.5 x 10^^^ erg 
cm~^ s~^, as is also the case for the directly selected hard 
X-ray sample. 

Another way to estimate the fraction of optically se- 
lected sources that are X-ray sources is to use a spectro- 
scopic data sample. Of the 554 / < 23 objects in the 4' 
radius region, 172 have measured redshifts or are spec- 
troscopically identified stars. The measured 2—10 keV 
fluxes of these galaxies are shown versus their absolute 
/-band magnitudes, M/, in Fig. |lo|b. Here K-corrections 
have been computed using the Coleman, Wu & Weedman 
(1980) spectral energy distribution (SED) for an early spi- 
ral galaxy. 



Of the 25 galaxies with -22.5 < Mj < -24, five have 
hard X-ray fluxes above 1.0 x 10"^^ erg cm~^ s~^. By 
contrast, none of the 29 galaxies with Mj fainter than 
—20 is an X-ray source at this level. If we correct for the 
incompleteness of the spectroscopically identified optical 
sample, we find that the fraction of optically luminous 
galaxies that are X-ray sources is O.OT^to'os- Thus, a very 
substantial fraction of the optically luminous galaxies are 
undergoing X-ray activity at any given time. 

6. SUBMILLIMETER AND RADIO PROPERTIES OF THE 
HARD X-RAY SAMPLE 

The new population of highly obscured, exceptionally 
luminous sources discovered by SCUBA (Small, Ivison, & 
Blain 1997; Hughes et al. 1998; Barger et al. 1998; Barger, 
Cowie, & Sanders 1999; Eales et al. 1999, 2000) appear 
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to be distant analogs of the local ultraluminous infrared 
galaxies (ULIGs; Sanders & Mirabel 1996). There is an 
ongoing debate on whether local ULIGs are dominantly 
powered by star formation or by heavily dust enshrouded 
AGN, and the same applies to the distant SCUBA sources 
(see, e.g., Trentham 2000). If the hard X-ray sources are 
highly absorbed, then the rest-frame soft X-ray through 
NIR radiation will be reprocessed by dust and gas, and the 
energy will appear in the FIR. At high redshifts (z ^ 1) 
the FIR radiation is shifted to the submillimeter. 

Barger et al. (1999) carried out a spectroscopic survey 
of a complete sample of submillimeter sources detected in 
a survey of massive lensing clusters. Only three of the 
17 sources could be reliably identified spectroscopically; of 
these, two showed AGN signatures. Thus, the possibil- 
ity that most SCUBA sources contain AGN remains open. 
Several authors (Almaini, Lawrence, & Boyle 1999; Gunn 
& Shanks 2000) have modelled the X-ray and submillime- 
ter backgrounds; they predict an AGN contribution to the 
SCUBA surveys at the level of 10 — 20 percent. 

The results of recent searches for submillimeter coun- 
terparts to Chandra X-ray sources have been mixed. In a 
study of two massive lensing clusters, A2390 and A1835, 
Fabian et al. (2000) identified three significant 2 — 7 
keV sources, but these were not significantly detected in 
the submillimeter (Small et al. 1998). Likewise, Horn- 
schemeier et al. (2000) did not see either of their 2 — 8 kcV 
sources in the ultradeep Hubble Deep Field SCUBA map 
of Hughes et al. (1998). In contrast, both of the 2-10 
keV sources detected by Bautz et al. (2000) in the A370 
lensed field are submillimeter sources (Small, Ivison, & 
Blain 1997). These two sources were previously identified 
spectroscopically as AGN (Ivison et al. 1998; Barger et al. 

1999) . The above mixed results probably reflect the fact 
that the 850 /zm flux limits obtainable with SCUBA are 
quite close to the expected fluxes from the obscured AGN, 
as we address below. 

The present study has the advantage that wide-area sub- 
millimeter (Barger, Cowie, & Sanders 1999 and the present 
paper) and extremely deep 20 cm data (Richards et al. 
2001) exist over the entire X-ray field of 57 arcmin^. An 
ultradeep 50 hr submillimeter map (Barger et al. 1998) 
also exists for one region of the field that contains two 
hard X-ray sources. 

6.1. Submillimeter Detection of a Hard X-ray Source 

With one exception, the hard X-ray sources are not de- 
tected in the submillimeter at the Scr level. The exception 
is source 8, an optically faint, highly absorbed (F = —0.06) 
X-ray source in the ultradeep submillimeter map. 

The bolometric FIR flux is related to the rest-frame 20 
cm flux through the well-established FIR-radio correlation 
(Condon 1992) of local starburst galaxies and radio quiet 
AGN. The SEDs of submillimeter sources are reasonably 
well approximated by the thermal black-body spectrum of 
the ULIG Arp 220 (Carilh & Yun 2000; Barger, Cowie, 
& Richards 2000), which is powered by star formation 
(Downes & Solomon 1998). We can therefore use the sub- 
millimeter to radio flux ratio to infer an approximate red- 
shift of 1.8 (with estimated redshift range 1.2 — 2.4) for 
source 8 (see Eqs. 2 and 4 of Barger, Cowie, & Richards 

2000) . The millimetric redshift estimation technique is ex- 
pected to hold for sources dominated by star formation or 



for radio quiet AGN. 

6.2. Spectral Energy Distributions of the Hard X-ray 

Sources 

The two quasars are the softest sources in the hard X-ray 
sample. In Fig. p4] a we compare their average rest-frame 
SED (filled squares) to the rest-frame SED of the z = 0.4 
radio quiet quasar PC 1543-1-489 (solid curve), normalized 
to approximately match the soft and hard X-ray data of 
the averaged quasars. 

Most of the sources in the sample show much harder X- 
ray spectra than the two quasars and are consistent with 
being obscured AGN. In Fig. [TT^ j we compare the rest- 
frame SED of source 8 (filled squares), based on the mil- 
limetric redshift of 1.8, with the SED of the local heavily 
dust obscured AGN NGC 6240 (solid curve), normalized 
to approximately match the hard and soft X-ray data of 
source 8. 

6.3. Bolometric FIR Luminosities Inferred from Radio 

Fluxes 

The sources with spectroscopic identiflcations are 
mostly at modest redshifts where radio observations pro- 
vide a more sensitive route for obtaining total FIR lu- 
minosities than do submillimeter observations (Barger, 
Cowie, & Richards 2000). Assuming that the sources 
are well described by the FIR-radio correlation (which 
holds for both star forming and radio-quiet AGN) given 
in Sanders & Mirabel (1996) with q ~ 2.35, then fpiR = 
8.4 X 10^** /20, where /20 is the flux at 20 cm. We calculate 
the FIR luminosity by relating it to the observed 20 cm 
flux, assuming a synchrotron spectrum with an index of 
-0.8; then 

Lfir = 4 tt 4 (8.4 X 10^^) /20 (1 + z)""'' 

where di, is the luminosity distance in cm and /20 has units 
erg cm~^ s~^ Hz^^. 

The bolometric FIR luminosities obtained by this 
method are given in column 5 of Table |^ and are plotted 
in Fig. 12 as open circles. Four sources are not detected 
above the 15 /iJy (3cr) limit in the 20 cm data and their lu- 
minosities are given as upper limits based on the 3a radio 
limit. For our sources without spectroscopic or millimetric 
redshifts, we nominally assume z = 2.0. It is interesting 
to note that three-quarters of the hard X-ray sources with 
radio detections in our sample would be classifled as LIGs 
(Lfib^ > 10" Lq; Sanders & Mirabel 1996) from their 
FIR luminosities, and half would be classifled as ULIGs 
[Lfir > lO^" Lq). 

6.4. Bolometric FIR Luminosities Inferred from 
Submillimeter Fluxes 

By assuming an ULIG SED, we can use the 850 /im 
data to estimate the FIR luminosity of source 8 and to 
place upper limits on the FIR luminosities of the remaining 
sources. For our non-quasar X-ray sources with redshifts, 
we estimate the FIR luminosities by scaling the NGC 6240 
luminosity by the ratio of the 850 /xm fluxes (after placing 
NGC 6240 at the redshift of the source). For NGC 6240 
we use the T = 42 K cold component determined by Klaas 
et al. (1997), which has Lfir = 2.0 x lO'^^ he5 erg 8"^ 
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Fig. 11. — (a) The average rest-frame SED (ufi,) of the two z = 2.56 quasars (filled squares) is compared with the rest-frame SED of the 
radio quiet quasar PG 1543+489 (solid curve). The radio to optical/NIR data on PG 1543+489 were taken from the compilation in PoUetta 
et al. (2000), the soft X-ray data from Laor et al. (1997), and the hard X-ray data from Vaughan et al. (1999). (b) The rest-frame SED of 
source 8 (filled squares) is shown in comparison with the SED (from Hasinger 2000) of the heavily dust obscured local AGN NGC 6240 (solid 
curve). 



in column 6 of Table 
squares. The radio in 



For our X-ray sources with no submillimeter detections, 
we use the 3cr upper limits on the submillimeter fluxes in 
calculating the Lpjn values. For our sources without spec- 
troscopic or millimetric redshifts, we nominally assume 
z = 2.0. Since above z ^ 1 the luminosity distance depen- 
dence and the K-correction approximately cancel, there 
is a simple translation of flux to bolometric luminosity 
independent of redshift; thus, our procedure should give 
accurate upper limits for the unidentified sources as long 
as they are not at low redshift (z < 1). 

For our two quasar X-ray sources, we estimate the FIR 
luminosities as above using the SED and FIR luminos- 
ity of PG 1543+489. For PG 1543-1-489 we calculate 
LpiR — 10.5 X 10'*^ /i65 erg s""'^ over the range 40 — 500 /im 
using the T = 46 K grey-body model given in PoUetta et 
al. (2000). 

The bolometric FIR luminosities or upper limits ob- 
tained from the submillimeter fluxes in this way are given 
3 and are plotted in Fig. O as filled 
erred FIR luminosities or upper lim- 
its are consistent with the submillimeter inferred FIR lumi- 
nosities or upper limits. It is clear from Fig. |l^ that more 
sensitive submillimeter observations would be required to 
detect even the higher redshifts objects. 

6.5. Bolometric Luminosity Ratios 

In Fig. ^ we plot the UV/optical to FIR luminos- 
ity ratios versus the logarithm of the soft to hard X-ray 
flux ratios. If the X-ray sources are hardened by line- 
of-sight opacities and the same material extinguishes the 
UV/optical light, then as the opacities decrease, we would 
expect both the UV/optical to FIR luminosity ratios and 
the soft to hard X-ray flux ratios to increase. For the 
spectroscopically identified sources there is a broad over- 
all trend in this sense but with large scatter. 

In Fig. Ijwe plot the ratios Lfir/ Lhx (open squares) 
and LoptT Lhx (filled diamonds) versus redshift. Only 
the two quasars and source 2 &t z — 1.320, with its cu- 
rious optical absorption line spectrum, are dominated by 
the UV/optical fight rather than by the FIR light. The 



two quasars are the softest sources in the sample and have 
X-ray photon indices (1.75 and 1.80) consistent with no 
absorption. Source 2 has a harder photon index (0.9) than 
the quasars; thus, the dominance of the UV/optical light 
here is perhaps less expected and may be due to substantial 
contributions of stellar light. For the remaining sources the 
FIR light dominates. 

6.6. Bolometric Correction 

If we assume that the FIR light in the hard X-ray sources 
is reprocessed AGN light, uncontaminated by star forma- 
tion in the galaxies, then we may use the present data 
to compute the bolometric correction from the hard X- 
ray luminosity to the bolometric light of the AGN. Ex- 
cluding the two quasars, which have an average (Lpip + 
Lqpt) / Lhx = 200, the remaining sources have an av- 
erage Loft/ Lhx — 19. If the sources with 3cr upper 
limits on L pip arc included at this level, then the average 
LpjpJ Lhx — 59, where we have used the radio inferred 
FIR luminosities, and the average {Lpjp + Lqpt)/ Lhx — 
78. If the sources with upper limits are instead assigned 
zero luminosity ratios, then the average Lpjuf Lhx = 56 
and the average {Lpm + Lqpt)/ Lhx = 74. 

The ratios are weighted to higher values by the two 
lowest redshift sources. If we instead weight the aver- 
ages by the hard X-ray fluxes, then the Lpip/ Lhx ratio 
drops to 33, the {Lpjp + Lqpt) /Lhx ratio to 42, and 
the Lqpt / Lhx ratio to 9. Because galaxy contamination 
of the optical light is a problem, the actual bolometric 
correction probably lies somewhere between the values of 
LpipJ Lhx and {Lpjp + Lqpt) / Lhx- As a conservative 
minimum, we adopt the Lpjp/ Lhx = 33 ratio for the 
bolometric correction in subsequent discussions while rec- 
ognizing that this ratio could be larger by a factor of two 
or more. 

The average Lqpt/ Lhx ratio for the hard X-ray 
sources is much smaller than the J5-band to 0.5 — 4 keV lu- 
minosity ratios seen in local early-type galaxies, whether 
we consider the soft components thought to arise from 
thermal emission in the gaseous halos (luminosity ratios 
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Fig. 12. — FIR luminosities versus redshift for the hard X-ray sample. Open circles were calculated from the radio data and the FIR-radio 
correlation. Filled squares were calculated from the submillimeter data and the NGC 6240 or PG 1543+489 SEDs and luminosities (see text 
for details). The spectroscopically unidentified sources, excluding source 8, are placed at z = 2; source 8 is placed at its millimetric redshift of 
1.8. Downward pointing arrows indicate 3cr limits for either the radio inferred FIR luminosities or the submillimeter inferred FIR, luminosities. 
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Fig. 13. — UV/optical to FIR luminosity ratios versus soft to hard X-ray flux ratios. Quasars are denoted by filled diamonds, AGN by 
filled triangles, 'normal' galaxies by filled squares, and spectroscopically unidentified sources by open squares. The latter, excluding source 8, 
are placed at z = 2; source 8 is placed at its millimetric redshift of 1.8. Upward pointing arrows indicate 3cr upper limits on the radio inferred 
FIR luminosities. 



in the range 150 — 40000) or the hard components which 
may be produced by X-ray binaries (luminosity ratios in 
the range 1200-8000) (Matsumoto et al. 1997). The much 
higher X-ray luminosities of the present sources relative to 
their stellar luminosities provide strong evidence that they 
are powered by AGN. 

It is also interesting to consider the total X-ray light. 
For an assumed photon index of F = 2, the total X-ray 
light is only weakly sensitive to the adopted energy range. 
For an energy range from 0.1 to 100 keV, the ratio of the 
total X-ray luminosity to the 2 — 10 keV luminosity is only 
a factor In(lOOO)/ ln(5) = 4.3. The total X-ray luminosi- 
ties of the sources, Lx, are still substantially smaller than 
the bolometric luminosities at other wavelengths, but an 
appreciable fraction of the AGN light (typically 1 — 20%) 
is emerging in the X-rays. Lbol = Lfir + Lqpt + Lx 



values are given in the last column of Table g; here again 
the radio determinations are used as the LpiR inputs. 

6.7. FIR and Hard X-ray Luminosity Comparisons 

By considering the ratio of the FIR luminosity, Lpm, 
to the hard X-ray luminosity, Lux, we can eliminate any 
redshift or cosmology dependence and make relative com- 
parisons of distant and local systems. Figure |l5| a shows 
the Lpin/ Lhx ratio versus Lhx for our data sample. 
Following Fabian et al. (2000), we also show on the fig- 
ure the data values for the local ultraluminous infrared 
galaxies Arp 220 and NGC 6240 (Klaas et al. 1997; Iwa- 
sawa 1999), the z = 0.4 radio quiet quasar PG 1543-1-489 
(PoUetta et al. 2000; Vaughan et al. 1999), and the radio 
loud quasar 3C 273 (Kim & Sanders 1998; Tiirler et al. 
1999). NGC 6240 hosts a powerful AGN and is highly ab- 
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Fig. 14. — Ratios of the bolometric optical (filled diamonds) and radio inferred FIR (open squares) luminosities to the hard X-ray lumi- 
nosities versus redshift. Spectroscopically unidentified sources, excluding source 8, are placed at 2: = 2; source 8 is placed at its millimetric 
redshift of 1.8. Downward pointing arrows indicate Scr limits for the radio inferred FIR luminosities. 



sorbed with an inferred column density of Nh ^ 2 x 10^"' 
cm"^ (Vignati et al. 1999). On the figure we have con- 
nected with a straight hne the values of Lhx calculated 
from the observed flux (Iwasawa 1999) and from the in- 
ferred intrinsic flux (Vignati et al. 1999) to illustrate the 
large effects that absorption can make for highly obscured 
AGN at low redshift. 

Using the submillimeter inferred FIR luminosity of 
source 8 and the radio inferred FIR luminosities of the 
other hard X-ray sources, we find that the values of 
Lfir/ Lhx for all the sources are comparable to that of 
NGC 6240. Also shown in Fig. |l5|a (open diamonds) are 
the luminous z > 4 radio quiet quasars observed in the 
submillimeter by McMahon et al. (1999) for which there 
are also X-ray detections (Kaspi et al. 2000; we converted 
the fi, at 2 keV values given in their Table 2 to 2 — 10 keV 
fluxes assuming F — 2). These quasars are in the region of 
3C 273 in Fig. HSa. The two gravitationally-lensed sources 
in the field of A^TO (open triangles) detected in hard X- 
rays by Bautz et al. (2000) and in the submillimeter by 
Small, Ivison, & Blain (1997) have higher Lpm/ Lhx ra- 
tios than the typical X-ray source. However, none of the 
hard X-ray selected sources remotely approaches the high 
Lfir/Lhx = 3.4 X 10^ value of Arp 220. 

Figure 15 b shows Lpjn/ Lhx versus Lhx for the above 
submillimeter sources and for submillimeter sources in the 
literature that have millimetric redshifts (Barger, Cowie, 
& Richards 2000 and § |^) and X-ray detections or lim- 
its (Hornschemeier et al. 2000 and § 0). If most of the 
submillimeter sources are star-formers like Arp 220, then 
much deeper hard X-ray observations would be required 
to detect them. 



7. A SUBMILLIMETER SELECTED SAMPLE: HARD X-RAY 
PROPERTIES AND CONTRIBUTION TO THE HARD XRB 

The foregoing section presented the submillimeter na- 
ture of the hard X-ray sources. We now invert the ap- 
proach and ask what are the hard X-ray properties of a 
submillimeter selected sample. We may address this by 
looking at the ensemble properties of submillimeter se- 
lected sources in the SSA13 field. There are twelve 3a 



source detections at 850 /im in the S2 and S3 chips within 
a 4.5' radius of the optical axis. The fluxes range from 
just over 2.3 mJy to 11.5 mJy. The positional accuracies 
for the submillimeter sources are relatively poor because 
of the large beam size. Six of the twelve submillimeter 
sources have 20 cm counterparts within 5", and the dis- 
persion of the offsets is 3.4" (see also Barger, Cowie, & 
Richards 2000). In order to allow for the positional un- 
certainty, we determined the 2 — 10 keV fluxes for each 
of the submillimeter sources in 10" diameter apertures. 
The results are shown in Fig. |l^. The only submillimeter 
source that is also a strong hard X-ray emitter is source 8 
from Table 1. The submillimeter source with the second 
strongest hard X-ray flux in Fig. |l6| is not in the hard X- 
ray sample but is a known soft X-ray emitter (source 28 
in Table 1 of MCBA). 

Because of the large apertures used, there is a signifl- 
cant probability of random overlap with a non-coincident 
X-ray source. In order to quantify this, we have again run 
Monte Carlo simulations using an equal number of sources 
at random positions. The la dispersion is 0.5 x 10~^^ erg 
cm~^ s~^ mJy^^. Only 5% of the simulations have values 
in excess of 1.2 x 10~^^ erg cm~^ s~^ mJy~^. In fact, the 
two sources that are contributing all of the X-ray signal 
have 20 cm radio counterparts that are closely coincident 
with the X-ray sources. Thus, both of these are real iden- 
tiflcations rather than chance contamination. 

The ratio of the total hard X-ray flux in the sample to 
the total submillimeter flux is 1.4 ± 0.5 x 10~^^ erg cm~^ 
s^^ mJy~^. We may use this ratio to estimate the fraction 
of the hard XRB that arises from submillimeter sources in 
this flux range. The EBL of the 2 — 10 mJy source popu- 
lation is 9.3 X lO'^ mJy deg~^ (Barger, Cowie, & Sanders 
1999), which would contribute 1.3 x 10~^^ erg cm~^ s""'^ 
deg-2 in the 2-10 keV band or 6% of the hard XRB. This 
contribution would rise to 19% to 27% if we normalized to 
the total submillimeter background of 3.1 x 10^ mJy deg~^ 
of Puget et al. (1996) or 4.4 x 10"* mJy deg^^ pixsen et 
al. (1998). 

However, nearly all of the X-ray signal from the sub- 
millimeter sample is coming from source 8. If this single 
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Fig. 15. — (a) T^firI ^HX versus Lux for the hard X-ray sample with rcdshifts. Quasars are denoted by filled diamonds, AGN by 
filled triangles, 'normal' galaxies by filled squares, and the optically faint source 8 by an open square. Also shown in the figure (see text for 
references) are the data values for the local ULIGs Arp 220 and NGC 6240, the radio quiet quasar PG 1543+489, and the brightest quasar in 
the sky, 3C 273. A straight line has been drawn between the absorption-corrected (higher X-ray luminosity) and observed NGC 6240 ratios to 
illustrate the large effects absorption can make at low redshift. The luminous z > 4 radio quiet quasars detected in both the submillimeter and 
X-ray are shown by the open diamonds. The two significant hard X-ray and submillimeter sources in A370 are shown by the open triangles, 
(b) Lpjjil Lux versus L^x for > 3cr submillimeter sources with spectroscopic or millimetric redshifts and hard X-ray detections or limits. 
Open circles denote the sources in the Hubble Deep Field and Flanking Fields, and open stars denote the sources in the SSA13 field (see text 
for references). Symbols for the comparison sources and for the sources significantly detected in hard X-rays have been filled in for clarity. 
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Fig. 16. - Hard X-ray versus submillimeter fiuxes for the submillimeter selected sample in the SSA13 X-ray field. The strongest hard X-ray 
emitter is source 8. The second strongest hard X-ray emitter is not in Table 1 but is a known soft X-ray emitter. 



source is removed, the total hard X-ray to total submil- 
limeter flux ratio drops to 0.6 ± 0.5 x 10""'^^ erg cm~^ s"""^ 
mJy^^. If we place the submillimeter sources at z = 2 
(consistent with the z — l — 'i spectroscopic and millimetric 
redshifts from, e.g., Barger et al. 1999 and Barger, Cowie, 
& Richards 2000), then the la lower limit on the ratio of 
LpiR to Lux would be approximately 1100. This lower 
limit is above the obscured AGN values and approaching 
that of Arp 220. It therefore appears that most of the sub- 
millimeter sources, at least above 2 mJy, are star formers 
with a small admixture of obscured AGN. Similar conclu- 
sions have recently been reached by Fabian et al. (2000), 
Hornschemeier et al. (2000), and Severgnini et al. (2000). 



8. A RADIO SELECTED SAMPLE: HARD X-RAY 
PROPERTIES AND CONTRIBUTION TO THE HARD XRB 

There are 107 20 cm sources above the 5cr radio limit of 
25 /iJy that lie on the S2 and S3 chips within a 4.5' radius 
of the optical axis. We have measured their 2 — 10 keV 
fluxes using 2.5" diameter apertures. Above a hard X-ray 
flux of 3 X 10""*^^ erg cm~^ s""'^, we recover only the 10 
hard X-ray sources with 5cr radio counterparts in Table 1. 
Above 1 X 10~^^ erg cm~^ s~^, 16 sources are detected. 
The contribution to the hard XRB from the entire radio 
selected sample is 5.9 x 10^^^ erg cm^^ s^^ deg^^ (26%), 
nearly all of which comes from the hard X-ray sources in 
Table 1. The small fraction (10 — 15%) of radio sources 
that are hard X-ray sources, and therefore AGN, is con- 
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sistent with expectations that the great bulk of the 20 cm 
sources at these faint fluxes are due to synchrotron emis- 
sion arising from supernova remnants in the interstellar 
medium of the galaxies rather than to nuclear activity. 
The radio morphologies, which will provide more insight 
into the nature of the hard X-ray sources, will be discussed 
in Richards et al. (2001). 

9. HARD X-RAY CONTRIBUTION TO THE EBL 

The summed contributions of our Chandra hard X-ra; 
sources to the EBL are shown versus wavelength in Fig 
as the filled diamonds. These are compared with measure 
ments of the EBL (solid curves) and with the integrated 
light from direct counts (large squares); the latter now lie 
close to the EBL at all wavelengths. 

For galaxies at z » 1 the X-ray power reprocessed by 
dust is expected to appear at submillimeter wavelengths. 
We computed the hard X-ray contribution to the sub- 
millimeter EBL by considering the average submillimeter 
properties of the whole X-ray sample and of selected sub- 
samples. Error- weighted sums of the submillimeter sources 
are given in Table together with the correspondingly 
weighted total hard X-ray fluxes for the sample and the 
ratios of the two. The total hard X-ray sample has a sig- 
nificant (> 3cr) submillimeter flux of 20 mJy. As expected, 
most of this arises in the z > 1.5 sources and in the uniden- 
tified sources. The total 850 /im flux of the sample cor- 
responds to an EBL of 4 x 10"^ mJy deg~^, or from 9 to 
13% of the 850 /im EBL, depending on whether the nor- 
malization of Fixsen et al. (1998) or Puget et al. (1996) is 
adopted. These values are similar to the obscured AGN 
model predictions of Almaini et al. (1999) and Gunn & 
Shanks (2000). However, it appears the the great bulk 
of the 850 /im EBL must arise in sources which are too 
faint at 2 — 10 keV to be seen at the current submillimeter 
threshold. 

We computed the hard X-ray contribution to the opti- 
cal/NIR EBL by summing the fluxes corresponding to the 
corrected 3" diameter aperture magnitudes at these wave- 
lengths. This light is dominated by the brighter galaxies, 
many of them with 'normal' galaxy spectra, so there is al- 
most certainly strong contamination of the estimate of the 
AGN contribution by the light of the host galaxies. The 
optical/NIR EBL should therefore be considered strictly 
as an upper bound. 

In summary, the hard X-ray sample contributes about 
10% of the light at both UV/optical and submillimeter 
wavelengths. Both wavelength regimes have similar to- 
tal bolometric light densities. However, because there is 
galactic light contamination in the optical/NIR, it is likely 
that the FIR is where most of the AGN light is emerging 
from the hard X-ray sources. 

10. BLACK HOLE MASS ACCRETION 

For nearby normal galaxies, there is an approximate em- 
pirical relation between the black hole mass, Mb/i, and the 
absolute magnitude of the bulge component of the host 
galaxy, M(bulge) (e.g., Magorrian et al. 1998; Wandel 
1999; Ferrarese & Merritt 2000; Gebhardt et al. 2000). 
Assuming that the relation also holds at high redshift, we 
can estimate the black hole masses of our 'normal' galaxies, 
subject to the uncertainties in translating our observed B- 



band magnitudes into Ms (bulge) that would worsen the 
empirical relation (Kormendy & Ho 2000). The Mb val- 
ues for our hard X-ray sample are plotted versus redshift in 
Fig. |l^. The inferred Mb/j values in our 'normal' galaxies 
are in the approximate range 5 x 10^ M0 to 10^ Mq. 

As discussed in § ||, luminous hard X-ray sources are 
common in bulge-dominated optically luminous galaxies 
with about 10% of the population showing activity at any 
given time. This preferential activity of luminous sources 
with massive bulges presumably reflects the relation be- 
tween the central massive black hole and the luminosity of 
the bulge. If the fraction of galaxies showing such behav- 
ior reflects the fraction of time that each galaxy spends 
accreting onto its massive black hole, then we require each 
such galaxy to be active for somewhere between 1 and 2 
Gyr. This is considerably longer than the theoretically es- 
timated accretion time of 0.01 Gyr for black hole fuelling 
by mergers (Kauffmann & Haenelt 2000) and therefore 
may suggest that the activity is being powered by smaller 
mergers or by internal flows within the galaxies. 

We may address the issue of the total accreted mass 
required to account for the luminosity in each galaxy by 
conservatively assuming that the hard X- ray flux- weighted 
average of Lpjji/ Lhx = 33 derived in § 6.6 is the ratio of 
the AGN's bolometric luminosity to its 2 — 10 keV luminos- 
ity; hereafter, we denote this bolometric correction by A. 
The mass inflow rate is eAMbhC^ / — Lbol — ALhx, 
where e is the efficiency for re-radiation of the accretion en- 
ergy. The total mass flow over the accretion time At ~ 1.5 
Gyr is 



AMbh = 10^ M, 



Our 'normal' galaxies have hard X-ray luminosities rang- 
ing from 0.1 to 50 x 10^^ erg s~^, for which the above 
accretion masses are 10^ to 5 x 10^ Mq. Thus, even for 
the maximum plausible efficiency, e ~ 0.1, we are seeing 
a substantial fraction of the growth of these supermassive 
black holes (Fabian & Iwasawa 1999). 

We may quantify this in a fairly model-independent 
way as follows. The bolometric surface brightness at the 
present time is simply related to the 2 — 10 keV XRB by 




Ibol = 0.85 A Ihx 



where the factor 0.85 is discussed in § 3.3. The bolometric 



EBL is related to the energy production as 



Ibol 



£ c-" pbh 

4 TT (1 -I- Z) 



where p^u is the universal mass density of supermassive 
black holes, and z is the mean redshift of the contributing 
sources. The (1 -|- z) factor reflects the adiabatic expan- 
sion loss. The hard X-ray flux-weighted mean redshift is 
z = 1.3 for the spectroscopically identified sources (includ- 
ing source 8) and z = 1.5 if we place the remaining sources 
at z = 2. 

We combine the two equations and normalize to the ob- 
served 2 — 10 keV background from Vecchi et al. (1999), 
Ihx = 7.6 x 10~* erg cm" 
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Fig. 17. — Filled diamonds show the contribution of our Chandra hard X-ray selected sample to the EBL at each of the measured wave- 
lengths. The error bars are smaller than the symbol size at all wavelengths. The dashed line shows the cosmic microwave background. 
Solid curves show the measured total cxtragalactic light in each wavelength regime: subiuillimeter from Puget et al. (1996) and Fixscn et 
al. (1998); optical/NIR from Bernstein et al. (1999) and Wright & Reese (2000): X-ray from Marshall et al. (1980), both at the original flux 
and re- normalized to the more recent higher flux determinations discussed in the text. Filled squares show the integrated contribution of the 
source counts: submillimeter counts to 0.5 mjy from Blaiii et al. (1999), but since these counts are quite uncertain, the well-established > 2 
mjy counts from Barger, Cowie, & Sanders (1999) are also shown {open square); optical/NIR from Williams et al. (1996) and Gardner et al. 
(1993); 1-2 keV X-ray from Table 2 of MCBA. 



Table 3 



Sample 


Total 


Total 


E/(2-10keV)/ 




/(850 urn) 


/(2 - 10 keV) 


^/(850 /xm) 




(mJy) 


(10-^® ergs s"^) 


(10"" ergs s"^ rajy"^) 



All 


20 ±6 


1800 


o.9j:? 


z < 1.5 


1±6 


840 


6.0?°2 


z > 1.5 (non-QSOs) 


7±3 


340 


o.5g:i 


QSOs 


5±2 


290 


o.6j:2 


Unidentified 


8±3 


600 


o.8j:i 




Pbh = 2x 10" 



where we have again conservatively used A ~ 33. Inclusion 
of the bolometric optical light would raise this estimate by 
about 40%. 

If we adopt a spheroid mass density of 10~^^ g cm~'^, 
which has a plausible multiplicative error of about 2 
(Cowie 1988: Fukugita et al. 1998), the black hole-bulge 
relation expressed in mass terms 

Mbh = 0.002 - 0.006 M(bulge) 

gives a black hole mass density of pbh = 2 — 6 x 10~^^ g 
cm~^. Thus, even for high-end estimates of the radiative 
efficiency (e ~ 0.1) we are seeing a significant fraction of 
the black hole accretion. Conversely, since we are proba- 
bly missing some fraction of the activity from sources that 
are too obscured to be seen in the 2—10 keV sample, we 
must have a high radiative accretion efficiency (e S> 0.01) 
or the required black hole mass density would be too high. 



We may go one step further and ask when the black 
holes formed, since the above equations also apply to the 
sources divided into redshift slices. We know directly from 
the spectroscopic obsc;rvations that at least 16% of the 
2 — 10 keV background arises below z = 1, 18% between 
z = 1 and 2, and 13% between z = 2 and 3. The remain- 
ing 21% of the 2 — 10 kcV light seen in the current sample 
most plausibly lies in the z = 1.5 — 3 range. Weighting by 
the mean redshifts in each interval, we find that at least 
10% of the observed black hole mass formation occurs at 
z < 1. This fraction could rise if there are more obscured 
hard X-ray sources missing from the 2 — 10 keV sample 
which preferentially lie at low redshift. 

11. CONCLUSIONS 

We carried out an extensive multi- wavelength observa- 
tional program to determine the nature of the hard X-ray 
background sources. Our major conclusions are as follows 

• In the 57 arcmin^ Chandra SSA13 field, we detected 
20 sources with 2-10 keV fluxes greater than 3.8 x 10~^^ 
erg cm~^ s~^. 

• We spectroscopically identified the 13 hard X-ray 
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Redshift 

Fig. 18. — Absolute S-band magnitudes versus redshift. 'Normal' galaxies are denoted by filled squares, AGN by filled triangles, and 
quasars by filled diamonds. The K-corrections have been computed using the Coleman, Wu, & Weedman 1980 SED for an early spiral galaxy. 



sources brighter than / = 23.5 mag; all arc in the red- 
shift range z = — 3. The spectra fall into three general 
categories: (i) 2 quasars, (ii) 5 AGN, and (iii) 6 optically 
'normal' galaxies. For the latter category, the AGN are 
either very weak or undetectable in the optical. 

• The soft to hard X-ray flux ratios of the hard sam- 
ple can be described by a rather narrow range of neutral 
hydrogen column densities from Nh = 2 x 10^^ cm~^ to 



3 X lO^'"' 



At most three of the 14 sources at z > 1 



have column densities above Nh = 3 x 10^^ cm~^, which 
suggests that we are seeing most of the obscured AGN in 
the present sample. 

• The hard X-ray sample is consistent with a constant 
surface density of ^ 400 sources per square degree per 
unit redshift. The redshift distribution is very similar to 
that of previous soft selected samples with similar limiting 
sensitivities. 

• The hard X-ray luminosities of the spectroscopically 
identified sources range from just over lO''^ erg s~^ to 
~ 10^^ erg s~^. All but two are more X-ray luminous 
than the local galaxy populations. 

• The integrated 2-10 kcV light of the 20 hard X-ray 
sources corresponds to an EBL of 1.34 x 10~^^ erg cm~^ 
s~^ deg~^, or between 58 to 84 percent of the hard XRB, 
depending on the XRB normalization. Most of this light 
is from the z <2 population. 

• The colors of the spectroscopically unidentified sources 
suggest that they are likely early galaxies at z = 1.5 — 3. 

• The 'normal' galaxies and AGN follow the upper en- 
velope of the star forming field galaxy population. 

• Between 4 and 12 percent of / < 23 galaxies with 
—22.5 > M] > —24 have hard X-ray fiuxes above 
1.0 X 10~^^ erg cm~^ s~^. Nearly all of these lie in the 
redshift range 2; = 0.3 — 1.5, for which the flux threshold 
corresponds to Lhx ^ 0.3 — 3 x 10^^ erg s~^. 

• Excluding the two quasars, the hard X-ray flux- 
weighted average bolometric optical to 2 — 10 keV luminos- 
ity ratio for the hard X-ray sources is low (9), providing 
strong evidence that the sources are powered by AGN. 



• Excluding the two quasars, the hard X-ray flux- 
weighted average bolometric correction from the 2 — 10 
keV luminosity to the bolometric light of the AGN is some- 
where between 33 and 42, depending on the optical light 
contamination by stars. 

• Of the 20 hard X-ray sources in our sample, only one 
(source 8) was significantly detected in the submillimeter. 
The millimetric redshift of source 8, obtained from the 
submillimeter to radio fiux ratio, is in the range 1.2 — 2.4. 
Its rest-frame SED is similar to that of the heavily dust 
obscured local AGN NGC 6240. 

• Bolometric FIR luminosities or upper limits inferred 
from the radio data using the FIR-radio correlation or 
from the submillimeter data using the NGC 6240 or 
PG 1543-1-489 SEDs and luminosities are consistent. 

• The > 2 mJy submillimeter source population con- 
tributes very little to the hard XRB (6%). 

• Excluding the one submillimeter source that had a 
significant hard X-ray detection, the FIR to hard X-ray 
luminosity ratio for the submillimeter selected sample has 
a IfT lower limit of 1100; this is above the obscured AGN 
values and approaching that of Arp 220. Thus, it appears 
likely that most of the > 2 mJy submillimeter sources are 
star formers. 

• The radio selected sample contributes 26% to the hard 
XRB, nearly all of which comes from the observed hard X- 
ray sources. 

• The hard X-ray sources contribute about 10% of the 
light at both UV/optical and submillimeter wavelengths. 
However, because there is stellar light contamination in 
the optical/NIR, it is likely that most of the AGN light is 
emerging in the FIR. 

• The masses of the black holes in our 'normal' galax- 
ies are estimated from the Magorrian relation to be in the 
range 5 x lO'^ M© to lO'' Mq. 

• Luminous hard X-ray sources are common in bulge 
dominated optically luminous galaxies with about 10% of 
the population showing activity at any given time. Since 
the hard X-ray emission is likely associated with episodes 
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of accretion onto the central massive black holes in these 
galaxies, the 10% may represent the duty cycle of galaxies 
that are active at any given time. Even with a high-end 
estimate of the radiative efficiency (e = 0.1), the black hole 
mass density required to account for the observed light is 
comparable to the local black hole mass density. 
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Fig. 1. — Chandra hard X-ray image of the SSA13 field. Tiic 20 significant hard X-ray source positions axe identified by the small circles. 
The large circle illustrates the 4.5' radius adopted for the analysis. 
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Fig. 2. — 9" X 9" _B-band thumbnail images of the hard X-ray sources in Table 1. The images arc from ultradeep data obtained with LRIS 
on Keck. The identification numbers arc as in Table 1, and the sources arc ordered from lower left to upper right by decreasing 2 — 10 keV 
flux. A circle of 1.5" radius, typical of the maximum positional uncertainty, is superimposed on each thumbnail. Redshifts, where available, 
are given in the lower left-hand corners of the thumbnails. North is up and East is to the left. For the more luminous objects we have shown 
the images at higher surface brightness to allow the positions of the objects with respect to the centers of the galaxies to be seen. 
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Fig. 3. — The Keck LRIS spectra of the 13 hard X-ray sources with secure redshift identifications. The spectra are shown in the rest-frame 
with source number and the measured redshift at the bottom of the frame. Dashed hues mark strong atmospheric absorption or emission 
features. 
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Table 1. 



# 




RA(2000) 


Dec(2000) 


/(2-10 keV) 
(lO-i*' cgs) 


/(0.5-2 keV) 
(10-1'' cgs) 


r 


B 


/ 


HK' 


AO 

(arcsec) 


2 


S(850/<m) 
(mJy) 


S(20 cm) 

(/"Jy) 


AO 

(arcsec) 


0(7) 


13 


12 


36.57 


42 


40 


2.9 


384 ± 34 


1480 ± 67 


1.54 


26.8 


24.0 


19.7 


(22.2) 


0.6 




1.1 ± 1.1 


21 




1 (35) 


13 


11 


59.36 


42 


39 


28.8 


273 ± 36 


612 ± 51 


1.01 


27.0 


23.0 


19.1 


[21.3) 


0.3 


1.048'' 


-3.9 ± 3.1 


76 


0.3 


2(9) 


13 


12 


35.68 


42 


41 


51.2 


232 ± 26 


408 ± 35 


0.90 


23.0 


20.8 


17.7 


(21.8) 


0.7 


1.320'' 


-3.9 ± 2.5 


5 




3 (24) 


13 


12 


15.25 


42 


39 


1.3 


190 ± 23 


986 ± 54 


1.75 


18.7 


17.9 


16.0 


(22.2) 


0.6 


2.5651 


3.5 ± 2.0 


74 


0.4 


4 (26) 


13 


12 


10.02 


42 


41 


29.3 


151 ± 20 


246 ± 28 


0.76 


22.5 


19.3 


16.6 


(21.9) 


0.1 


0.212" 


2.1 ± 1.9 


14 




5 (29) 


13 


12 


6.58 


42 


41 


41.8 


138 ± 20 


96 ± 18 


0.05 


23.8 


20.5 


17.8 


(21.2) 


0.2 


0.696^ 


-0.3 ± 1.8 


44 


0.3 


6 (19) 


13 


12 


22.39 


42 


38 


14.4 


116 ± 20 


622 ± 45 


1.80 


21.6 


19.8 


17.8 


(22.5) 


0.6 


2.565^ 


2.2 ± 1.6 


29 


0.5 


7 (10) 


13 


12 


34.44 


42 


43 


9.0 


78 ± 15 


3 ± 13 


< -2.5 


21.9 


18.8 


16.1 


(21.2) 


0.9 


0.241" 


1.0 ± 1.8 


80 


0.6 


8 (15) 


13 


12 


28.27 


42 


44 


54.6 


57 ± 14 


33 ± 14 


-0.06 


25.7 


24.3 


21.4 


(22.3) 


0.6 




2.0 ± 0.6 


39 


0.4 


9(5) 


13 


12 


37.12 


42 


43 


21.3 


50 ± 12 


106 ± 22 


1.03 


26.9 


24.7 


20.5 


(21.7) 






0.6 ± 1.8 


22 




10 (16) 


13 


12 


26.08 


42 


37 


36.7 


50 ± 15 


170 ± 25 


1.43 


25.0 


22.7 


19.5 


(21.7) 


0.7 


1.427" 


0.0 ± 1.7 


18 




11 (22) 


13 


12 


19.98 


42 


42 


21.5 


49 ± 12 


196 ± 25 


1.53 


24.5 


22.3 


19.6 


(22.1) 


0.3 


2.415" 


1.3 it 2.0 


25 




12 (12) 


13 


12 


31.27 


42 


39 


1.6 


48 ± 13 


20 ± 9 


-0.40 


23.1 


20.2 


17.7 


(22.0) 


0.3 


0.585" 


0.0 ± 1.4 


18 




13 (32) 


13 


12 


4.23 


42 


41 


13.8 


48 ± 14 


26 ± 10 


-0.17 


> 27.6 


> 25.9 


> 21.9 


(21.9) 






-1.6 ± 1.7 


9 




14 (6) 


13 


12 


36.71 


42 


38 


44.0 


46 ± 14 


39 ± 12 


0.28 


27.0 


25.5 


20.7 


(21.5) 






1.0 ± 1.5 


31 


0.7 


15 (1) 


13 


12 


40.25 


42 


39 


35.7 


43 ± 14 


79 ± 17 


0.93 


24.7 


23.3 


20.9 


(21.9) 


0.6 


2.625" 


1.3 ± 1.2 


-11 




16 (18) 


13 


12 


22.48 


42 


44 


50.9 


42 ± 12 


151 ± 26 


1.49 


26.3 


> 25.9 


20.8 


(22.6) 






0.6 ± 0.7 


17 




17 (8) 


13 


12 


36.04 


42 


40 


44.2 


42 ± 12 


23 ± 9 


-0.06 


25.5 


23.8 


20.6 


(22.1) 


0.2 




0.6 ± 1.2 


77 


0.2 


18 (3) 


13 


12 


39.36 


42 


42 


49.1 


39 ± 11 


205 ± 28 


1.82 


18.9 


15.9 


13.5 


(21.5) 


0.4 


0.110" 


0.7 ±2.2 


91 


0.2 


19 (30) 


13 


12 


6.45 


42 


41 


26.6 


38 ± 12 


8 ± 6 


-1.72 


19.8 


16.7 


14.2 


(21.5) 


0.5 


0.180" 


1.5 ± 1.8 


91 


0.2 
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